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ABSTRACT 



We present a direct detection of the stellar halo of the edge-on SO galaxy NGC 3957, using ultra-deep VLT/VIMOS V and R images. This is 
achieved with a sky subtraction strategy based on infrared techniques. These observations allow us to reach unprecedented high signal-to-noise 
ratios up to 15 kpc away from the galaxy center, rendering photon-noise negligible. The Icr detection limits are R = 30.6 mag/arcsec" and 
y = 31.4 mag/arcsec^. We conduct a thorough analysis of the possible sources of systematic errors that could affect the data: flat-fielding, 
differences in CCD responses, scaling of the sky background, the extended halo itself, and PSF wings. We conclude that the V - R colour of 
the NGC 3957 halo, calculated between 5 and 8 kpc above the disc plane where the systematic errors are modest, is consistent with an old and 
preferentially metal-poor normal stellar population, like that revealed in nearby galaxy halos from studies of their resolved stellar content. We 
do not find support for the extremely red colours found in earlier studies of diffuse halo emission, which we suggest might have been due to 
residual systematic errors. 
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1. Introduction 

Galaxy stellar halos contain fundamental clues about galaxy 
assembly. ACDM models predict steUar halos to be ubiquitous 
and dominated by old metal-poor populations, characterized by 
significant subst ructures in the form of tidal debris from ac- 
creted satelhtes jBullock & JohnstonlEool lAbadi et alJl2006t 
iFont et al. I l2006h . Conversely, in classical dissipative collapse 
models, halos are ex pected to exhibit sm ooth age and metal- 
licity gradients (e.g., lEggen et al. Due to the extreme 

faintness of these envelopes, the value of galactic stellar halos 
as key tests of galaxy formation theories has however not yet 
been fully realized. 

Since the middle of the nineties, our view of galaxy 
haloes has changed dramatically with the discovery of sub- 
structures arotind the Milky Way (s -p- 



20081 iMcConnachie et al.1 120091: iTanaka et alj I2OI0I) . Stellar 
streams were also detected around NGC 5907 and NGC 4013 
jMartfnez-Del^ado et aDl2008lE00 9). It is nevertheless not yet 
totally clear how ubiquitous halos are and how often mergers 
intervene in their building-up. 

Resolving individual stars is mainly limited to a small sam- 
ple of galaxies within or close to the Local Group, with a 
few exceptions. iMouhcine et al. I (l2005allbl) resolved the stars 
near the red giant branch tip in the outskirts of eight nearby 
galaxies with the HSTAVFPC2, although it remains unclear 
for some of their fields w hether they sampled the g alaxy outer 



discs or halos. Simila r ly. | Mpuhcine et al.l ( 120071) : llbata et al 
(l2009l) : lReikubaet al.l (l2009h analysed the extra-planar stel- 



Ibata et al. 
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Maiewski et al. 1999; Newbergetal. 2002; Gilmore et al. 
2002: Rocha-Pinto et al. 2003: Vivas & Zinn 2006.: .Grillmair 



2006[|Belokurov et alj2007l:ljuric et alj2008l:lBeIl et alj2008l) 

Similar signatures of tidal destruction of satellites b y their mas- 
sive hosts were found in the halo in M31 (e.g. 



Ibata et al 



200ll:lFerguson et alJl2002Hlbata et ai]l2007l:lRichardson et al 



* Based on observations obtained at the ESO Very Large Telescope 
(VLT) in the Program 077.B0145(A) 



lar popul ations in NGC 89 1 u sing three HST/ACS pointings. 
Recently ' Barker etaL (l2009h conducted the first wide-field 
ground-based survey of the red giant branch population in the 
outskirts of M81 using Suprime-Cam on the 8-m Subaru tele- 
scope. They detected a faint, extended structural component be- 
yond the galaxy's bright optical disc sharing some similarities 
with the Milky Way's halo, but also exhibiting some important 
differences. 

Searches for extended low surface brightness diffuse halo 
emission were carried out around several more distant external 



galaxies. The first detection was announced by ISackett et al 
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(11994 aroun d the edge-on Sc ga laxy NGC 5907 and late r 
confirmed bv lLequeux et al.l d 19961) and lLequeux et alJ (Il998l) . 
These two subsequent works additionally reported an unex- 
pected reddening of the halo colours, increasing with the dis- 
tance to the galaxy center These results were contested by 
Zheng et alJ (119991) . who noticed that artifacts in the star- 



maskin g procedure could b e responsible for the apparent red- 
dening. IZibetti et al. I (12004 used a stacking technique to com- 
bine the images of 1047 edge-on galaxies selected from the 
Sloan Digital Sky Survey. They detected the presence of a mean 
luminous halo, whose shape was clearly rounder than the disc. 
Their data suggested a correlation between the halo and galaxy 
luminosities. Surprisingly, their halo colours were redder 



than th e reddest known elliptical galaxies. IZibetti & Ferguson 
( 12004 reached similar photometric depths in a study of a single 
edge-on galaxy in the Hubble Ultra Deep Field, where they also 
found anomalously red colours, not reproduced by any conven- 
tional stellar populations model and getting redder with larger 
radii. |de^ Jong (2008) however discussed how scattered galaxy 
light from extended point spread function tails could have been 
underestimated in these studies and thus lead to spurious detec- 
tions. 

The controversy over the existence and nature of stellar ha- 
los around external galaxies reflects that these types of obser- 
vations are extremely challenging from a technical standpoint. 
Errors in flat-fielding and sky subtraction can easily mask or 
contort real signal. In order to overcome these limitations and 
to understand galactic halo properties, we designed a new strat- 
egy designed for galactic halo surface brightness profiles. It is 
based on two principles: i) we target a galaxy, which has an 
angular size that is much larger than the PSF extent and ii) we 
carry out very accurate sky subtraction, with a chopping tech- 
nique similar to that used in near-IR imaging. 

This paper is organized as follows: Section |2] presents our 
observational technique. Section[3]describes the data reduction 
procedure leading to ultra-deep sky-subtracted V and R images 
of NGC 3957. Section]?] presents the radial surface brightness 
profiles derived for the halo of NGC 3957. Section]5]discusses 
the most important possible sources of errors, focusing on the 
systematics. Section ]6] interprets the V - R colour profile. Our 
conclusions are summarized in Sect. ]7] 



2. Observational strategy 

We selected the early-type (SO) galaxy NGC 3957 
(RA(J2000):11 54 01.5 ; DEC (J2000) : -19 34 08 ; 
v/,c;=1637 km/s). NGC 3957 is seen nearly edge-on, its 
exact inclina t ion an gle (88 + 1 degrees) was measured by 
Pohlen et all (]2004 . Its apparent size (3.1' x 0.7') is large 
enough to allow us to integrate flux over large regions when 
determining the surface brightness profiles. At the same time, 
it is small enough to fall entirely in a single VIMOS CCD 
chip. Taking Ho = 73 km/s/Mpc, at the distance of NGC 3957 
(22.42 Mpc), 1 arcmin=6.52 kpc. Another important point 
is that the suiToundings of NGC 3957, which were used to 
determine the background level and shape, are devoid of any 
bright objects. 



Table 1. Journal of the observations. Half of the exposures of 
each night were taken with the galaxy located in CCDl and the 
others with the galaxy in CCD2 (see Fig.]2]l. 



Date 
Gjjj/mm/dd) 



Exposure Filter 



2006/04/24 
2006/04/25 
2006/04/27 
2006/04/28 
2006/04/29 
2006/05/25 



24 X 190s 
24 X 190s 
12 X 190s 
24 X 190s 
24 X 190s 
24 X 190s 



2006/05/27 
2006/05/28 
2006/05/29 
2006/05/31 



12 X 450s 
12 X 450s 
12 X 450s 
14 X 450s 



V 



The choice of VIMOS was motivated by the opportunity to 
have a wide field of view distributed over four different CCDs, 
of 2048 X 2440 pixels each. With a resolution of 0.205" per 
pixel, each CCD covers 7' x 8'. They are separated by a gap of 
2' (see Fig.]2]i. Unfortunately, the choice of a guide star satisfy- 
ing VIMOS specifications had the consequence that the guide 
probe fell on two of the four CCDs. In all data, it covers more 
than 50% of either CCD3 or CCD4 (see CCD4 in Fig.]2]l, leav- 
ing those two chips unusable for our analysis. 

We gathered R and V-band images, with total exposures 
of 22500s and 25080s (i.e., ~ 6 and 7 hours), respectively. 
Exposures were split into sequences of 190 s (/?-band) and 
450 s (V-band). The journal of the observations is presented 
in TablelH 

The design of VIMOS allowed us to place NGC 3957 in 
one CCD while integrating the sky on the other one. This strat- 
egy offered two major advantages: (/) On one hand, the back- 
ground level was estimated in a region located far away from 
the galaxy, thus avoiding possible contamination by faint, ex- 
tended galaxy components, (ii) On the other hand, a 2-D image 
of the sky was obtained under the exact same observing condi- 
tions as those of the galaxy. 

Exposures were dithered following the patterns shown in 
Fig. JT] This allowed us not only to build master skies from 
empty field frames, nicely removing any bright object, but also 
to smear out any inhomogeneities in the background, increas- 
ing the flatness of the combined master skies. 

As already mentioned, we used two different pointings. We 
will use the following terminology throughout: Pointing 1 cor- 
responds to a set of exposures placing the galaxy in CCDl 
and dithered following the patterns of FiglT] Pointing 2 has 
the same dithering pattern as Pointing 1 , but the galaxy is now 
placed in CCD2. The exposures on the empty field in CCDl 
(when the galaxy is in CCD2) and CCD2 (when the galaxy 
is in CCDl) are used to determine the background level and 
shape. 
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(150,150) 




Fig. 1. Left: dithering pattern for a set of 12 consecutive /?-band frames. The position (0,0) indicates the initial position of NGC 
3957 in one CCD. Right: dithering pattern for a set of six consecutive V-band exposures. Units are given in pixels. One pixel is 
0.205". 
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Fig. 2. Geometry of the VIMOS field of view, composed of four CCDs, prior to sky subtraction. The galaxy was alternately put 
in CCDl and in CCD2. The black feature in CCDS and CCD4 is due to the guide star camera. 
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Fig. 3. Example of the master sky in R (left panel) and V (right panel) The cuts are +3cr, the large scale variations are < 1 %. 



3. Data reduction and sky subtraction 

We only describe at length the steps differing from the clas- 
sical methodology for the reduction of images. These steps 
are critical to reach a surface brightness level down to ~ 30 
mag/arcsec^. 

Bias and flat-field corrections were done in a classical man- 
ner. Five featureless biases were taken each night. Their mean 
ADU/pixel was subtracted from all frames of the same night. 
When possible, a series of five twilight flat fields were taken. 
For each of these nights, we composed a master flat, corre- 
sponding to the average of the five flats, with a 2-cr clipping 
boundary, cr being calculated with the gain and readout noise of 
the CCDs. Since flat-fields were not taken each night, the same 
master flat was used for several subsequent nights. We will dis- 
cuss the possible uncertainties arising from the flat-fielding in 
Sect.O 

The most crucial step of the data reduction is the sky sub- 
traction. One master sky per night and per CCD was created 
as the average (plus 2-cr clipping boundary) of all dithered 
frames which did not contain the galaxy. For the 7? -band, 12 
sky images per night and CCD were available. Sigma clip- 
ping of the dithered images proved sufficient to remove any 
bright star in the field. For the V-band, only six frames were 
available. Consequently, we had to mask the bright stars be- 
fore dithering to remove them completely from the master sky 
frames. Figure [3] provides an example of a master sky for the 
two photometric bands. The cuts are chosen to be +2% of the 
mean background level, which is three times the standard de- 
viation. The large scale structures represent <1% of the mean 
background level. They are of the order of the uncertainties in 
flat-fielding, which are discussed in Sect. 15.11 

For the sake of clarity, we will now use the following 
nomenclature (see Fig.|4|: 

- MSI: Master sky of CCDl, taken from the empty field 
frames of Pointing 2, 

- SKYli: individual CCDl image of Pointing 2, which does 
not contain the galaxy and from which MS 1 was created. 



- GALli! individual CCDl image of Pointing 1, which con- 
tains the galaxy, 

- MS2: Master sky of CCD2 with Pointing 1. Analogous to 
MSI. 

- SKY2i: individual CCD2 image of the empty field with 
Pointing 1. Analogous to SKYli, 

- GAL2i: CCD2 individual galaxy image of CCD2 with 
Pointing 2. Analogous to GALli, 

The index i goes from 1 to 12 for the 7? -band or from 1 to 6 
for the V-band. 

We could not directly use the sky frames obtained at the 
same time as the galaxy. Indeed, galaxy and background would 
then have been observed with different CCDs. Small differ- 
ences in the flat-fielding between CCDl and CCD2 could in- 
troduce artificial intensity gradients in the final images. We 
needed to subtract skies of Pointing 2 from the galaxy frames 
of Pointing 1 (and vice versa) obtained during the same night. 
In practice, for each night we subtracted MSI from GALli 
and MS2 from GAL2i, i.e. the master skies were observed at a 
slightly different time than the galaxy, but with the same CCD. 
Although full dithering patterns for both pointings were com- 
pleted in about one hour, the sky level was not necessarily sta- 
ble in this time interval. Actually, we measured a variation in 
the mean background level of 10-15% from the first to the last 
exposure for a given pointing sequence. Two coiTections were 
therefore necessary: 

(1) For each GALli (GAL2i), MS1(MS2) had to be scaled 
to the level of the sky at the time at which the galaxy was ob- 
served. We had direct access to this information via the contem- 
poraneous exposures SKY2i (SKYli). To determine the scaling 
factor of each GALli, MSI was multiphed by the ratio 



< SKYli > 
<MSl> ' 



(1) 



where <MS1> is the mean intensity, in ADU/pixel, of a 
large central region (x; 2M pixels, i.e., 22 kpc-radius region 
centered at ~50kpc from the galaxy center) of the master sky 
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frame and <SKY2i > is the mean intensity in the same region 
of the sky field frame, observed simultaneously with GALli. 
Our surface-brightness limit is ~ 31.4 mag/arcsec^ in V and 
~30.6 mag/arcsec- in R. These limiting magnitudes correspond 
to the integrated sky noise in 1 arcsec^, or equivalently, to the 
fluxes leading to a signal-to-noise rat io of 1 (28 .2 mag/arcsec- 
in R and 29.0 arcsec^ in V at S/N=3). ^in et al. (2005) report 
a V-band halo surface brightness in M31 of 31-32 mag/arcsec^ 
between ~ 40 and 55 kpc along the minor axis. Since this area 
of M3 1 is contaminated by tidal debris streams, the actual sur- 
face brightness of the underlying halo is somewhat lower (e.g. 
Ibata et alJ (120071) V Assuming that the halo of M3 1 is identical 
to that of NGC 3957, it corresponds to the level of noise in our 
sky frames. (2)Since < SKY2i > and < MS 1 > are measured 
on two distinct CCDs, we had to correct /3 for their different re- 
sponses. It was done by calculating the ratios of the intensities 
of the bias-subtracted, non-normalised flatfields of CCDl and 
CCD2. Namely, each /3i had to be multiplied by a sensitivity 
correction 



CCD 1 



CCD 2 



< FFl > 
< FF2 > ' 



(2) 



for the scaling of the MS 1 frames and by a ' for the scaling 
ofMS2. 

Errors in a have opposite effects in CCDl and CCD2 re- 
spectively, leading to a systematic additive offset in the results 
when considering each of the two CCDs separately. This can 
be used to fine-tune a removing any systematic difference be- 
tween CCDl andCCD2. 

The complete sky subtraction is finally described by 



FINli = GALli -MSI - a -/in 



and 



A2 



FINli = GALli -MSI-—, 

a 



(3) 



(4) 



where FINli and FIN2i are the final sky subtracted frames 
of CCDs 1 and 2. A schematic view of the scaling method is 
shown in Fig.|4] 

The 132 (50) /?-band (V-band) background-subtracted im- 
ages were finally aligned and combined, and the final images 
are calibrated in the Johnson-Cousins photometric system. One 
CCD frame per night was calibrated in V and R with the stan- 
dard stars observed during the same period. The final V and R 
combined images were then scaled to the level of their respec- 
tive reference V and R frames. Hence, the final frames were 
multiplied by the factor flux(single)/flux(combined), where 
flux(single) and flux(combined) were measured from several 
bright stars in the reference and combined frames, respectively. 
These photometrically calibrated images were then corrected 
for extinction using Ar = 0.124 and Ay - 0.153 following 
Schlegel et al.l (Il998h . The extinction varies by 0.006 mag rms 
within a radius of 6 arcmin around NGC 3957. Considering 
random errors alone, the 1 -cr final accuracy of our calibration 
is 0.022 mag, in R, and 0.036 mag, in V. 



galaxy : 


rames 




empty fields 








a 




GALl. 
1 




SKY2 

1 1 



pointing 1 



master sky 

MSI from 
SKYl. 



galaxy frames 



GAL2. 



Fig. 4. Explanation of the scaling scheme, prior to the image 
coaddition. The master sky of CCDl (MSI) is taken as a ref- 
erence. The / index refers to the individual exposures within 
a sequence. The factors fSi scale MSI to match the sky level 
of each empty fields of CCD2 (SKY2i), observed on the same 
night. The factor a corrects for the difference in sensitivity be- 
tween the two CCDs. It is determined from the intensity ratios 
of the flatfields of CCDl and CCD2. 




Fig. 5. Cone shaped regions used to integrate the light on the 
final coadded V and R images. North is to the top, east to the 
right. 
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4. Surface brightness profiles 



In order to trace the galaxy signal down to very faint levels 
with a sufficiently high signal-to-noise ratio, we integrated flux 
in 90° wedges, running along the minor axis of the galaxy 
(Fig. lU. Prior to this binning, the foreground stars were care- 
fully masked, as the background galaxies. The radial width of 
each bin was chosen as a trade-off between the spatial res- 
olution of the galaxy profile and the desired signal-to-noise 
per bin. The narrowest bins have an annular width of 2 pix- 
els (0.04 kpc) and are located in the central parts of the galaxy, 
while the widest ones have a radial extent of 100 pixels (1 kpc). 
The resulting surface brightness profiles in V and /? -bands, ob- 
tained by taking the average of the flux in all unmasked pixel 
in each radial ring, are shown in Fig. |6] where the dashed lines 
indicate the detection limits per square arcsecond. The sp atial 
sampling is five times finer than that of IZibetti et al. (2004) in 
the inner galactic regions, while in the outer parts, the spa- 
tial sampling of th e two studies is the same. However, unlike 
Zibetti et al.l [ (12004 . who get a signal-to-noise of 1 for the out- 
ermost bins at ~12-14 kpc, we are not limited by photon noise 
in these parts (S/N~20). Systematic errors (see Sect. |5]l start 
dominating over photon noise at r=8-10 kpc, where our S/N is 
still w 40 in y and R . 

We estimated the photon noise cr,^, per bin by adding in 
quadrature the photon noise of all N pixels in a cone-shape bin 



(5) 



where the sum goes over all non-masked pixels. The noise 
per pixel, /, on the sky-subtracted and coadded frame is 



(Ti = 



— +0- 



sky^ 



(6) 



where Ft is the flux in pixel /, g is the electron-to-adu con- 
version factor, and cr,^,, is the standard deviation of the pixels 
in the sky background. This calculation gives typical error bars 
of 0.05 mag in regions located at about 10 kpc from the centre. 
The error bars remain smaller than 0. 1 mag from 10 kpc to the 
edges of the frames at 15 kpc. These small eiTor bars are too 
small to be seen on Fig.|6] An exception is the west-side of the 
V-profile. This is the faintest of the profiles and reaches the de- 
tection limit already at ~12 kpc. Therefore, the outermost data 
points beyond this distance are not considered further. 

5. Possible systematic errors 

Whilst our detection of the NGC 3957 stellar halo is not lim- 
ited by photon noise, it may still be affected by a number of 
systematic errors that we now review and try to quantify. This 
is of particular importance when trying to interpret the radial 
changes in the colour profiles. 

5.1. Flat fields 

Errors in flat-fielding may cause artificial gradients in the sky 
background and, in turn, affect the actual halo colour gradi- 
ent. Since our skyflats were not taken each night, a possible 
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Fig. 6. V- (black) and /?-band (red) surface brightness pro- 
files along the minor-axis of NGC 3957. Errors due to photon 
noise are typically <0.05 mag within 10 kpc and <0.1 mag for 
r > 10 kpc, hence too small to be plotted. The systematic er- 
rors are not displayed for clarity and can be seen in Fig. [8] The 
dashed lines indicate the Icr detection limits in the two bands. 



source of error is a the temporal variability of the flat-fields 
over a period of several nights. In order to estimate this vari- 
ability, we compare our normalized /?-band flat-field obtained 
at the beginning of the observing run (24 April 2006) with the 
one corresponding to the end of the observing run (25 May 
2006). The changes during this one month period do not ex- 
ceed 1.5% across the whole field of view. Since we used the 
same flat-fields for the sky and the galaxy frames (see Fig.|4|l, 
eiTors in the flat-fields propagate linearly into eiTors in the final 
profiles. Consequently, a 1.5% deviation in the flat-field results 
in a 1 .5% deviation in the the surface brightness. This translates 
into a maximum change of shape of 0.01 mag, i.e., completely 
negligible compared with the photon noise. 

5.2. Extended PSF tails 



In a recent work, de Jond ( 2008') discusses how the extended 
wings of the point spread function (PSF) can significantly con- 
taminate the measurement of galaxy halo light if not properly 
accounted for. 

By targeting a very low redshift galaxy with a size extend- 
ing over several arcminutes, we ensured that even the most 
extended parts of the PSF did not affect the galaxy's sur- 
face brightness profiles at all. Indeed, the PSF in our coadded 
frames has a FWHM~1", which is more than 100 times smaller 
than the exten t of the measured galaxy halo. For comparison, 
de Jongl ( 2008 ) deals with PSFs that are about ten times smaller 
than the observed galaxy. 
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Fig. 7. Mean flux level as measured in a 100 x 100 pixel empty 
box at the boundary of the individual, sky-subtracted frames. 
The measurements are for the /?-band. Black circles are for 
CCDl, while the red circles are for CCD2. There is no sys- 
tematic difference between CCDs, indicating that a is properly 
determined. 



In addition, we checked that the wings of the PSF have a 
limited size. At 1.8 arcsec (corresponding to 0.2 kpc at the dis- 
tance of NGC 3957) the flux in our PSF drops to 0.72% of its 
central intensity in R, and to 2.55% of the central intensity in 
V. The integrated flux of bright standard stars through apertures 
of growing radius showed that there was no measurable flux in 
the PSF wings akeady at 15" away from the core. Finally, we 
scaled the PSF to the central intensity of NGC 3957 and mea- 
sured the flux in its wings. At 50 pixels, i.e., 1.1 kpc from the 
PSF centre, the level of the flux was below O.OIADU, i.e., not 
measurable. It is clear that the wings of the PSF have no effect 
on our surface brightness profiles. 

5.3. Difference in tfie CCD responses: a -factors 

The a coefficient described in Sect. [3] (Eq. |2]l accounts for the 
difference in sensitivity between the two CCDs. The determi- 
nation of this factor for each single epoch was carried out in 
two steps. 

1- We initialized a first guess value by dividing CCDl and 
CCD2 flatfields by pairs. 

2- We applied a small correction to this initial value by 
noticing that the scaling of the sky frames subtracted from 
CCDl is proportional to a, while the scaling for CCD2 is pro- 
portional to a"'. As a consequence, a slight error in a is imme- 
diately seen as a systematic flux offset between the two CCDs, 
which is minimal when the value of a is such that the sky sub- 
tracts equally well in CCDl and in CCD2. In order to evaluate 
this offset and to correct for it, we determined in each individ- 



ual sky-subtracted frame the mean flux level in a 100 x 100 
pixels empty region close to the image boundary. Our final 
measurements are shown in Fig. [7] after correction; there is no 
residual offset between the mean a factors of CCD 1 and CCD2 

The amplitudes of the corrections applied to a were of 
the order 1% at most. This translates into changes of about 
0.01 mag/arcsec^ in the central regions of the galaxy and to 
0. 1 mag/arcsec^ beyond 1 3 kpc. The maximum change in slope 
of the profile, due to an erroneous estimate of a, is therefore at 
most 0.1 mag/arcsec^ in the faintest and most external regions 
of the galaxy field. 

5.4. Sky-level scaling: /3 -factors 

The determination of the mean sky level in each individual im- 
age is affected by shot noise, which in turn leads to an over- or 
under-subtraction through erroneous /3 factors (Eq. [T]). While 
quantifying this error for each image taken separately is im- 
possible, we used the full data set to estimate the amplitude of 
the fluctuations with time. This is equivalent to measuring the 
scatter of the /3 factors for a fixed a. 

The amplitude of the random fluctuations of the scale fac- 
tor j6 can be estimated from the scatter of the points in Fig. |2l 
The standard deviation of the 132 7? -band images is 4.66 ADU, 
while the uncertainty on the mean /3 is 0.41 ADU. Similarly, the 
standard deviation of the 50 V-band frames is 3.75 ADU, while 
the uncertainty of the mean value of /3 is 0.53 ADU. These val- 
ues provide good estimates of an upper limit on the errors on /? 
in V and R. 

A systematic error on /3 may have important consequences 
on the colour profile of the galaxy, in particular if the amplitude 
of the error is different in V and in R. In order to estimate this 
effect, we artificially introduced an over- or under-subtraction 
of the sky by 2cr in both bands. This corresponds to a shift of 
2 X 0.41 ^ 0.8 ADU in Randlx 0.53 -1.1 ADU in V. 

It translates into a variation of ~0.05 mag at 25.5 
mag. arcsec^, ~0.2 mag at 27 mag. arcsec^ and ~0.5 mag at 28 
mag. arcsec^ in the 7?-band. Similarly, the shifts in the V-band 
are ~0.05 mag, ~0.4 mag, and more than ~1.0 mag at 26, 28, 
and 30 mag. arcsec^, respectively. 

5.5. Extended halo 

The last possible source of error is the presence of a very ex- 
tended halo signal in the SKYlj and SKY2i frames used to 
scale MS2 and MS 1 before subtraction (see Fig.|4]i. This rescal- 
ing indeed assumes that the CCD frames, which do not contain 
NGC 3957, are far enough away from the galaxy to be free 
from any residual halo light. 

Since the reduction and scaling procedures are the same in 
the V and R filters, the systematics causing a positive (nega- 
tive) shift in V would also cause a positive (negat ive) shift in 



R. Fro m their radial minor axis profiles of M31, llrwin et al 



(120051) derive a power law surface brightness profile following 



I(r) oc r , beyond 20 kpc. We used this relation to evalu- 
ate to what extent SKYli and SKY2i might be contaminated 
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H (kpc) H (kpc) 

Fig. 8. The two-component fits to the minor axis profiles in the V and R -bands. The observed profiles are shown with solid dots 
with attached error bars (see Sect. |5]l. The last point of the V-West profile is not shown, neither is it used in the fits, because it 
is obviously dominated by errors (see Sect.|4]i. The red lines indicate the de Vaucouleurs law, which best fits the bulge of NGC 
3957. The red solid line indicates the region of the fits within 5 kpc, while the dashed line shows its extrapolation to larger radii. 
The blue lines show the power-law component of the fit at radii larger than 5 kpc. 



by extended halo light from one CCD on the other Our sky 
frames are located at a minimum of twice the distance from 
the galaxy cen t re to t he edge of our galaxy profile. Assuming 
the Irwin et al. I (l2005h M3 1 halo surface brightness can be con- 
sidered representative, the halo light in the CCDs we used to 
measure the sky should be at least five times fainter than at the 
edge of the CCD that contains the galaxy. 

This translates into an error of 1 ADU/pixel, which is a very 
conservative value, corresponding to about 40% of the mea- 
sured /?-band flux and 60% of the flux in the V-band at 13 kpc 
from the galaxy centre. This modifies the V and R East-side 
magnitudes by ~0.03 mag at 5 kpc, ~0.2 mag at 10 kpc, and 
~0. 3-0.4 mag (R) respectively ~0.5 mag (V) at 15 kpc. 



5.6. Total error budget 

To summarize, the main possible sources of errors that can af- 
fect the V and R surface brightness profiles of NGC 3957 pre- 
sented here are i) the uncertainties related to the scaling of the 
sky-level through the y6-factors and ii) the possible overestima- 
tion of the actual sky background level, caused by contamina- 
tion by a genuine extended stellar halo signal in the regions 
used to build the sky frames. This latter factor can only act as 
sky over-subtraction and therefore will contribute only to the 
upper limit of the error bars. Our surface brightness profiles 
are presented in Fig. |8] along with their total error bars. 



6. Properties of the stellar halo 

6. 1. Detection in V andR 

As shown in Fig.|6] we clearly detect light up to about 15 kpc 
away from the centre of NGC 3957. In order to evaluate the 
structure of this luminous component and, more importantly, to 
see whether one can resolve it into more than a single structure, 
we performed ID-fits to the V and R profiles. We considered 
only the photon noise error bars, since they are the only statis- 
tically meaningful quantities to be taken into account inx^ pro- 
cedures. We used the IDL routine mpfit, which solves the least- 
squares problem with the Levenberg-Marquardt technique, to 
conduct the fits. The inner 0.5kpc region was not considered 
due to the presence of the dust lane. We first started with a de 
Vaucouleur law form, chosen to best represent the bulge of this 
SO galaxy. Figure |8] shows that while a r'^^ function with an 
effective radius of 0.29 kpc provides an excellent description 
of the inner profile in both in V and R of NGC 3957, it consid- 
erably under-predicts the light beyond 4-5 kpc. Beyond 5 kpc, 
the excess light is well-represented by a power-law with index 
-2.76 + 0.43. A pure exponential fit with a scale length of a few 
kpc always provides a worse fits, pa rticularly in the j^-b and, for 
which it was essentially discarded. Ilbata et al. I (l2009h applied 
a de Vaucouleur law with an effective radius of 1.55 kpc to 
NGC 89rs extra-planar light, but this did not provide a good 
description of the light detected here. 



Pohlen et al.l (120041) present 3D thick/thin disc decomposi- 



tions for a sample of eight SO galaxies, including NGC 3957. 
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R (kpc) 



Fig. 9. The R-hand vertical profile of NGC 3957 extracted at 
a position 5 kpc along the maj or axis. The observ ations are 
displayed with open circles. The Pohlen et al. ( 20041) fits of the 
thin and thick discs are indicated with a dashed red line, while 
our pure exponential fits are shown with a plain blue line. 



They derive an /?-band vertical scale height for the thick disc 
component of 2.3 kpc, which, by their definition, corresponds 
to a normal exponential vertical scale height of 1.65 kpc. In or- 
der to compare this value with our own, we needed to extract 
a vertical profile some distance from the minor axis in order to 
eliminate the influence of the bulge. Figure |9] shows the V and 
R West and East profiles of NGC 3957 extracted at a position 5 
kpc along the major axis. Fitting the thick disc between 2 and 
4 kpc, we found an exponential scale height of 1.27 kpc, which 
agrees well with that found by Pohlen et al.l ( 20041) . consider- 
ing the different ways in which the two studies have extracted 
profiles and their different photometri c depths and s patial res- 



olutions. This fit and also that from Pohlen et al. I (l2004 



overlaid in the figure. Figure|9]also underscores the need for an 
additional component beyond a thick disc at vertical distances 
greater than ~4 kpc. 

6.2. Colour profile 

Figure [TO]displays the minor-axis colour profile of NGC 3957. 
The apparent strong reddening at large galactocentric distances 
must be analysed in light of the possible systematic errors pre- 
sented in Sect. |5] As seen in Fig. |6] the NGC 3957 surface 
brightness profile is traced down to about 1.2 mag above the 
detection limit both in V and R bands. 

We now analyse the west side of the galaxy (negative val- 
ues of the radii) as an illustration of the different elements at 
play. The total error bars in V and R are small and similar from 
the galaxy centre, up to about 8.5 kpc (A < 0.05 mag). At 
9.5 kpc, the error in V is 0.2 mag larger than the one in R. 
This difference keeps increasing with the radius, reaching 0.95 
mag at 10.5 kpc. This effect is due to the strong decrease in 
V-flux rendering the sky subtraction very critical, i.e., a given 
error on the sky subtraction has much larger impact on the 



> 





R (kpc) 

Fig. 10. TheV-R colour profile of NGC 3957. The edges of 
the vertical bars show the maximum possible systematic errors, 
which dominate over the photon noise, as estimated from the 
uncertainties in the sky-level scaling /? factors (Sect. 15.4b and 
considering possible extended halo (Sect. 15.51 ) 



V-band fainter profile than on the R-hand one. To be quan- 
titative, at 8.5 kpc, the V and R fluxes are at the level of ~4 
ADU/pixel and ~4.5 ADU/pixel, respectively. At 9.5 kpc, we 
measure only ~2 ADU/pixel and ~3.5 ADU/pixel in V and R, 
respectively. Figure [TO] clearly shows that the larger the error 
bars due to systematics, the steeper the rise inV - R. This sug- 
gests that this reddening may not be a genuine property of the 
stellar/dust properties of NGC 3957 halo, but rather an artifact 
due to the sharp drop in the V flux. The sky subtraction be- 
comes very uncertain and leads to biased results when the halo 
surface brightness approaches the detection limit by less than 
two magnitudes, as indicated in Fig.|6] This is the case at a dis- 
tance of 10 kpc from the galaxy centre in V, and at 15 kpc in 
the /?-band. 

The west-side mean colour profile has V - R ~0.50 +0.02 
from 2kpc from the centre of the galaxy, up to 5 kpc and then 
marginally rises to 0.70 +0.2, at 8 kpc. As discussed earlier, 
the error bars are not independent from each other and rather 
reflect systematic errors. They all act in the same direction, ei- 
ther increasing or decreasing the overall V - R profile. The up- 
per boundary (i.e., connecting the upper limits of the error bars) 
corresponds to an error of lADU/pix on the true sky level, or 
equivalently to a 3cr error (1.2 and 1.6 ADU in R and V, re- 
spectively) on the scaling of the sky. Although very modest, 
these values have noticeable consequences on the shape of the 
profiles, due to the faint luminosity levels considered here. In 
regions where the V - R gradient is steeper than 0.1-0.2 mag 
per kpc, the systematic errors dominate. The analysis of the 
east side of NGC 3957 follows the same line, albeit the mean 
colour is redder by ~0.05 mag with respect to the west-side, 
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even in the regions where the systematic errors are negligible. 
The small colour asymmetry between the east and west sides of 
NGC 3957 is intriguing and could possibly reflect fluctuations 
due to substructures within the galaxy's halo. 

We now investigate whether the inner halo colour of NGC 
3957 is indeed compatible with those of normal stellar popula- 
tions, and conduct comparisons with a few galaxies for which 
the halo stellar population could be resolved and analysed. We 
decomposed the halo population into the same number of sin- 
gle stellar populations (SSPs) as there are bins in the published 
metallicity distribution function (MDF). These SSPs were gen- 
erated with metallicit ies equal t o that of the bin mean using 



the Y isochrones of iKim et al.l (l2002h . that extend down to 



[Fe/H]=-3.8, wifli [a/Fe]=0.3, a Salpeter IMF and a 13 Gyr 
age. When necessary the isochrones were interpolated to rep- 
resent the metallicities of the MDFs. The resultant SSPs were 
then summed in luminosity, weighted by the fraction of light 
contributed by each MDF bin. 

We first derived the expected integrated V - R colour 
of the Milky Way halo from the metallicity distribution of 



Carollo et al.l (12009 ). They find that from a vertical distance of 
5 kpc up to 9kpc above the Galactic plane, which corresponds 
to the regions we sampled in the present study, the halo stellar 
metallicities span [Fe/H]=-3.0 to -1 with a small fraction of 
stars reaching [Fe/H]=-0.5. The distributions at 5.5 kpc and 8.5 
kpc, peaking at [Fe/H]=-1.5 at and [Fe/H]=-2, respectively , 
lead to y - R= 0.49 and 0.47. This narrow range of colours is 
due to the progressive insensitivity of the isochrones to changes 
in metallicity at low [Fe/H] . 

As to M31, the complex web of stellar streams prevented 
us from deriving a priori a unique colour for its halo. We in- 
stead examined the possible range of values based on earlier 
works providing complete metallicity distributions. The metal- 
licity distribution derived by Ibata et al. (2007) from the diffuse 
U ght of M31 halo tran slates into V - ~ 0.51. The M2 field 
of lDurrell et al.l (1200 ih . located at 20 kpc along the M3 1 minor 

axis, peaks at [M/H] 0.5 with a long tail of more metal-poor 

population down to -2.5, leading to an i ntegrated V -R ~ 0.54. 



In contrast, the spectroscopic study of iKoch et al.l (l2008i) re 



veals a strong metallicity gradient with a peak at [Fe/H] 1 . 

for galactocentric distances between 20kpc and 40kpc and at 
[Fe/H]~ -2 beyond 40kpc, with V-R -0.47. 

Now turni ng to slightly more distant galaxies 



Mouhcine et al.l (12007!) analyse the 1.5 first magnitude of 
flie RGB stars in the halo NGC 891, at -9.5 kpc from the 
galactic plane, centred on V - / ~1.5 and derive a peak 
metallicity of [Fe/H]~ -1, i.e, close to an integrated colour 
o f y - ~ 0.5. Som e of the galaxies in the sample of 



Mouhcine et al.l (12005 ah have metallicity distribution peaks at 

[M/H] 0.6. Their full metallicity distributions results in an 

integrated V -R ~ 0.54. 

In summary, V-R from ~0.45 to ~0.6 constitutes the range 
of expected integrated colours of old and preferentially metal- 
poor galactic halos. The halo of NGC 3957, with V-R from 
0.5 to 0.7 mag, between 5 and 8 kpc from the galaxy centre, 
where the stellar halo component is clearly dominating and the 
systematic errors are still modest, agrees fairly well with these 
numbers. Beyond these vertical distances, the colour redden- 



ing coiTelates with increasing errors. The order magnitude of 
the former being similar or smaller than that of the latter sug- 
gests that the colour gradient is due to the uncertainties in sky 
subtraction at very faint flux levels. 



7. Conclusions 

We obtained ultra-deep optical VLT/VIMOS images of 
NGC 3957, a nearby edge-on SO galaxy. The total exposure 
time of six hours in R and 7 hours in V allowed us to reach the 
limiting magnitudes of R = 30.6 mag/arcsec^ and V - 3 1 .4 
mag/arcsec^ in the Vega system. 

We devised a new observational strategy based on infrared 
techniques, which takes advantage of the large field of view 
of VIMOS. By observing NGC 3957 alternatively in two dif- 
ferent CCDs, we were able to create high signal-to-noise and 
stable sky frames, which were used to carry out accurate sky 
subtraction, without any assumption on the spatial shape and 
flux level of the sky background. The observations allowed us 
to reach unprecedented high signal-to-noise ratios up 15 kpc 
away from the galaxy centre, rendering photon-noise negligi- 
ble. They enabled the clear detection of the stellar halo of NGC 
3957 at distances above 5 kpc. 

We conducted a thorough analysis of the possible sources 
of systematic error that could affect the vertical surface bright- 
ness and colour profiles: flat-fielding, differences in CCD re- 
sponses, scaling of the sky background, extended halo, PSF 
wings. While a r'^"^ function with an effective radius of 0.29 
kpc provides an excellent description of the inner V and 7?-band 
minor axis profiles of NGC 3957, it under-predicts the light be- 
yond 4-5 kpc. Above 5 kpc, the galaxy light profile requires an 
additional power-law component with index -2.76 + 0.43. 

The most secure part of the NGC 3957 halo colour profile 
falls in the range V -R ~ 0.5 - 0.7 mag. This colour is compat- 
ible with the properties of nearby galaxy halos as revealed by 
the investigation of their resolved stellar populations and with 
that of the Milky Way. An apparent strong reddening is seen in 
the outer parts (r > 8 kpc) of the V-R profile, similarly to ear- 
lier works on galaxy halo surface brightness profiles. However, 
our analysis indicates that this reddening may not be a genuine 
property of the stellar/dust properties of NGC 3957 halo, and 
simply reflects the impact of systematic eiTors in sky subtrac- 
tion at extremely faint flux levels. It is possible that previously 
published reports of extremely red colours in external galaxy 
halos could have resulted from similar effects. 

Future imaging programmes benefiting from larger fields of 
view (~ 0.5 - 1 .0 degree on one side) will allow the systematics 
in halo emission searches to be controlled even further. Existing 
facilities to caiTy out the experiment include the SuprimeCam 
on the Subaru telescope and the VST (VLT survey telescope). 
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